Superfluid phases of ultracold Fermi gases on a checkerboard superlattice 
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We analyze the ground-state phase diagram of two-component Fermi gases loaded into a two-dimensional 
checkerboard superlattice, i.e. a double-well optical lattice, potential within the BCS mean-field theory. We 
show that, by coupling the two s-wave sublattice superfluid order parameters, a checkerboard potential gives 
rise to a Hamiltonian that has the form of a two-band superfluidity with three (two intraband and an interband) 
nonlocal order parameters. We study the evolution of these order parameters as a function of particle filling, 
interaction strength and checkerboard potential, and find that the system always prefers the 0-phase solutions, 
i.e. the phase difference between sublattice order parameters is 0, but never the 7r-phase one. In addition, we 
find that the ground-state of the system undergo a superfluid-normal quantum phase transition at half fillings 
beyond a critical checkerboard potential C, the threshold of which is precisely determined by the magnitude 
of the order parameter at C — 0, and that the normal state rapidly turns into a checkerboard insulator as C 
increases. 

PACS numbers: 05.30.Fk, 03.75.Ss, 03.75.Hh 



I. INTRODUCTION 



Optical superlattices UJj-|4j] have been of interest to the cold 
atom community for a long time now, and they have recently 
regained more interest in the past few years since they may 
allow for the studies of topological phases of matter in the 
cold-atom context 1 5-8] . One way to realize a checkerboard 
super optical lattice is to superimpose two independent optical 
standing waves differing in period by a factor of two, e.g. A/2 
and A, and tunable intensities and relative phases. In particu- 
lar, a double-well optical lattice can be produced by arranging 
the shorter-wavelength (A/2) lattice potential in such a way 
to split the each potential well of the longer- wavelength (A) 
lattice into two. In addition, the energy difference between 
the wells of the resultant double-well potential can also be 
controlled by tuning the relative phase of the optical poten- 
tials GUI. 

In spite of serious challenges in producing fermion su- 
perfluids in earlier optical lattice experiments [9-13], there 
is some recent experimental evidence for superfluid, metal- 
lic and insulating phases B14l - tl6ll (see also the recent re- 
view 111711 ). Motivated by these experiments, in this paper, we 
investigate the ground-state phases of two-component Fermi 
gases loaded into a two-dimensional checkerboard superlat- 
tice. For this purpose, we study a Fermi-Hubbard type lat- 
tice model which includes, in addition to the usual nearest- 
neighbor hopping and onsite (attractive) interaction, an onsite 
energy difference between sublattice sites, i.e. a staggered 
checkerboard potential. We note that the the phase diagram 
of the Bose-Hubbard versions of such a model have recently 
been studied for the hardcore [181 and softcore [191 bosons. 



Our main findings, within the single-band tight-binding 
BCS mean-field theory, are as follows. First, we show that 
the s-wave sublattice order parameters, which are momentum 
independent in the original Hamiltonian, are coupled by the 
presence of a checkerboard potential, and this gives rise to a 
Hamiltonian that has the form of a two-band superfluidity with 
three (two intraband and an interband) nonlocal (momentum- 
dependent higher partial waves) order parameters in the basis 



where the single-particle Hamiltonian is diagonal. We study 
the evolution of these order parameters as a function of par- 
ticle filling, interaction strength and checkerboard potential, 
and found that the system always prefers the 0-phase solu- 
tions, i.e. the phase difference between sublattice order pa- 
rameters is 0, but never the 7r-phase one. In addition, we 
show that the ground-state of the system undergo a superfluid- 
normal quantum phase transition at half fillings beyond a crit- 
ical checkerboard potential C, the threshold of which is pre- 
cisely determined by the magnitude of the order parameter at 
C = 0, and that the normal state rapidly turns into a checker- 
board insulator as C increases. 

The rest of the paper is organized as follows. We introduce 
the single- and many-body Hamiltonians in Sec. [Ill an d derive 
the single-quasiparticle/hole excitation spectra of the system 
as well as the complete set of self-consistency (superfluid or- 
der parameters and total and imbalance number) equations. 
We solve the resultant equations in Sec. [TTT1 and give a de- 
tailed analysis of the obtained results mentioned above. A 
brief summary of our main findings is given in Sec. |IV] 



II. HAMILTONIAN 



It is well-established that Hubbard-type discrete lattice 
models can be used to capture the physics of cold atoms 
loaded into optical lattice potentials J20H25I1 . For exam- 
ple, much of the theoretical predictions based on the sim- 
plest Bose-Hubbard model J26I1 have been successfully ver- 
ified with ultracold Bose gases loaded into optical lattices. 
The prime examples are the realizations of superfluid and 
Mott insulator phases as well as the transition between the 
two II20h25I1 . Motivated by this success, here we study a 
Fermi-Hubbard type lattice model to analyze the physics of 
ultracold Fermi gases loaded into a checkerboard superlattice 
potential, as described next. 



A. Single-particle problem 

Let us first discuss the single-particle problem on a 
two-dimensional checkerboard superlattice which consists 
of two interpenetrating square sublattices as illustrated in 
Fig. [TJ Within the tight-binding approximation, the hop- 
ping Hamiltonian for such a lattice can be written as Hq = 

-EieojeAo-WjS.aCfaa^i where i a >P} = ( A > B ) la " 
bels the sublattices and a — (f , J,) labels the two components 
of the Fermi gas, ti a ,jp :(T is the tunneling (hopping) ampli- 
tude of a fermions between lattice sites ia and jf3, and cj 
(ciaa ) operator creates (annihilates) a a fermion at lattice site 
ia. In this paper, we limit ourselves to the simplest model 
where ti a ,jp,a — ti<, is nonzero for all nearest-neighbor lattice 
sites but otherwise. However, we also allow the possibility 
of alternating hopping elements with a checkerboard pattern 
where £ may take two values. In particular, let's assume ti a 
(t 2a ) corresponds to the hopping amplitude from the sublat- 
tice A to the sublattice B in the positive (negative) x and y 
directions. We note that the next-nearest-neighbor hoppings 
are suppressed by an order of magnitude compared to t a in a 
typical atomic setting IU7I1 . 



the dj directions (see Fig. [TJ where d = y/2a, 
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FIG. 1: (Color online) The checkerboard superlattice consisting of 
interpenetrating A and B square sublattices is sketched. Here, a is 
the lattice spacing, and (e x = ax, e y = ay) and (di = e x 
e y , d2 = e x + e y ) are the primitive unit vectors used to construct 
the single-particle problem. 

The eigenvalues of such real-space hopping Hamiltonians 
can be obtained by taking advantage of the discrete transla- 
tional symmetry, and transforming them to the Fourier (mo- 
mentum) space. Therefore, we assume a single-band descrip- 
tion, and introduce c iaa = (1/-/&Q Z)k Q eiBz c k„ao-e lk °' ri , 
where Ma = Mb = M/2 is the number of sublattice sites, 
k Q is the momentum, 1BZ is the corresponding 1st Bril- 
louin zone, and r; is the position of the lattice site i. Using 
this transformation, we rewrite the Hamiltonian in momentum 
space as H Q = £ ko . (e kcT cl Al7 c kBa . + H.c. ) , where 



£ko 



2 (ti,r + t 2a ) cos(fcid/2) cos(k 2 d/2) 
2i{t 1(7 - t 2<T ) cos(kid/2) sm(k 2 d/2) (1) 



operator creates (annihilates) a a fermion with momentum 
k = (k x ,k y ) on a sublattice, and H.c. is the Hermitian 
conjugate. Thus, the single-particle dispersion relations are 
£kor = H e ko-|, where r = (+, — ) labels the two bands. Note 
that, since the translational symmetry is doubled in real space, 
the original k space is halved but the number of bands is dou- 
bled, in such a way that this dispersion recovers the usual re- 
sult, i.e. Eke- = — 2t a [cos(k x a) + cos(k y a)], in the absence of 
a sublattice structure when t\ a — t 2a = t a . 

In addition to the hopping part, we include an onsite 
checkerboard lattice potential, i.e. an alternating energy 
off-set between sublattices, which is given by H c \> = 
- c Y,ie a ,a7ac\ a<7 c laa , where C > is its strength and 
7 Q = +1 (—1) for a = A (B). This term lowers (raises) 
the onsite energy of the A (B) sublattice sites by C. Us- 
ing the Fourier transformation described above, we rewrite 
the total Hamiltonian He = Hq + H c \, in k space and di- 
agonalize it, le ading to H c = Eka^kmtL^. where 
£kcra = 7a\/|ekcr| 2 + C 2 are the single-quasiparticle/hole 
dispersion relations in the presence of a checkerboard poten- 
tial, and bL (b kaa ) operator is the new quasiparticle cre- 
ation (annihilation) operator in the transformed basis. Note 
that £k CT Q reduces to e kcrr when C — as expected. We 
note that the Ck QCT and b kaa operators are related via a Bo- 
goliubov transformation c k Ao = u ka Ab k Aa + UkcrsbkBa and 
Cks CT = v^AbkAa+VkaBbkBa, where u kaa anduko-a are the 
components of the eigenvector that corresponds to the eigen- 
value £kera ■ The eigenvectors are orthonormal in such a way 
that Mkcra/wktra = £kcr/(C + e kaa ). Having discussed the 
single-particle Hamiltonian, next we move on to the many- 
particle problem. 



B. Many-particle problem 

For the many-particle problem, the effects of local (onsite) 
and attractive interparticle density-density interactions can be 
taken into account within the BCS mean-field approximation, 
which is known to work well for weak interactions at all tem- 
peratures and even for moderate interactions at T = 11271 — 
29]. For this purpose, we introduce sublattice-dependent su- 
perfluid order parameters A^ and Ab with s-wave symme- 
try, as defined by A a = -(2g/M)J2 ieot (c ia ^c ia f) in real 
or -(2g/M)Y /k (c-kaiCkat) in k s P a ce, where g > is 
the strength of the interaction and (• • • ) is a thermal aver- 
age. Therefore, the interaction contribution to the k-space 
Hamiltonian can be written as H m f — (M/g) J2 a l^a | 2 /2 + 
J2 ka (A* c_kaiCkat + H.c.) . We also introduce a spin- 
dependent chemical potential /i CT term to the Hamiltonian, i.e. 

H n = - E*g q ,<t Av4kx<W in real or - J^kaa Av4 a(J c kQ0 . 
in k space, which allows us to fix the number of a fermions 
independently of each other. 

Thus, the total many-body mean-field Hamiltonian H = 



H + H c b 



H 



ml 



Hfj, for the checkerboard super- 



is in general a complex number. Here, fcj = k • di/d cor- 
responds to the projections of the momentum vector along 



lattice can be compactly written in k space as, H = 

(M/ ff )EJA Q | 2 /2~2E k ^i + EkV'kDkV'k,whereVk = 



( c kAt' c icBt' c -ks;, c_ k AiJ denotes the fermionic operators 
collectively, and the Hamiltonian matrix Dk is 



Dl 



i t -C 


e kt A A 


f* 

fc kf 


-Hf + C A B 





A s H~ C - £ -H 


A* 


-e^k, H + C 



(2) 



In this paper, we consider equal hoppings for j" and \. 
fermions, i.e. ti-f = t^ — tg, leading to ekf = e -ki = £ k- 
The single-quasiparticle/hole excitation spectra £k s of the in- 
teracting system are given by the eigenvalues of this Hamilto- 
nian matrix, and they can be compactly expressed as 



E h 



-h + ft 




(3) 



X = 4|e k | 2 (|A A | 2 + |A B | 2 - 2|A A A S | cos$) + 16e k /x 2 
+ (\A A \ 2 - \A B \ 2 ) (|A A | 2 - |A S | 2 + 8fiC) , (4) 



where s = (1,2,3,4) with ft = ft = -ft = -ft = 1 
labels the eigenvalues, h — (/if — /xjJ/2 is the difference 
and fi = (/z-j- + H\)/2 is the average chemical potential, 
e k = -\/|ek| 2 + C" 2 , and $ = 4>a — 0b is the phase difference. 
Here, we assume A a — |A Q |e^°. Note that the particle-hole 
symmetry of the Hamiltonian implies simultaneous transfor- 
mation of A — s- B and /i — > — /i. Note also that, after setting 
A^ = Ab = Ao in the C = limit, this expression recovers 
the usual result E^ s = —h + ft\/(|ek| — /i) 2 + |Ao| 2 , which 
is doubly degenerate since the original k space is halved. 

Using the quasiparticle/hole excitation spectra, we obtain 
the corresponding mean-field thermodynamic potential £1 for 
the total Hamiltonian H as 



y a ks 



1 + tanh 



(%H 



JE^ s -2^ Mj 



(5) 



written as 

|Aa| 



.9 

|A 5 | 

3 



^^tanb/^ 1 " 



2M <^ 9|A^ 

ks 



— E 



2M 



ks 



2T 

-^ks 



9S ks . , 

■— - — - tann , 
d\A B \ \ 2T 



2M ^ 

ks 



dEk s ( Ek s 

— - — tanh — — - 

d/j, \ 2T 



S^ dEk s ( E\^s 



2M ^ dh 

ks 



2T 



(6) 

(7) 
(8) 
(9) 



Some of these partial derivatives are long and not par- 
ticularly illuminating, and therefore none of them are 
shown. Similar to the expressions above, the remaining 
(fifth) phase-difference equation can be written as = 
£ ks (<9£ks/,9$)tanh[£ks/(2T)], a nd since dX/d<$> = 
8|ek| 2 |A,iA.B| sin$, we immediately conclude that $ is ei- 
ther or 7r. However, our numerical results suggest that 
the 7r-phase solution is never realized for the particular 
model Hamiltonian that we consider in this paper. Note 
that these k-space summations over the 1BZ can be con- 
verted into the k-space integrations via X)k/(^ 1 '^ 2 ) — 
[M/(8^ 2 )] /: w /: w f(x/d, y/d)dxdy. 

Before attacking this problem via numerical means, we 
would like to gain some physical insight. For this purpose, 
we rewrite the mean-field Hamiltonian H m f in the basis of 
b^Aa and &kBo- operators, i.e. the one where the single- 
particle checkerboard Hamiltonian is diagonal. Up to a con- 
stant term, the resultant four terms can be compactly written 
as Zlk Q /3(- F ka / 36_ka|6k/3t + H - c -)> where the coefficients are 
given by F^ afi = A A u-^ ia u^p + A* B v- k i a v ktl3 . It is easy 
to show that the intraband coefficients F^ a ^p vanish when 
C — >• as one may expect. This immediately reveals that 
H m f has the form of a two-band superfluidity Hamiltonian 
as long as C ^ 0. Note that the coefficients F^aa, F^bb 
and F-kAB correspond, respectively, to the k-dependent A in- 
traband, B intraband and AB interband superfluid order pa- 
rameters, and all of them are spin singlet with even parity, i.e. 
-Fka/3 = F-ka/3, as one may expect. Therefore, the starting s- 
wave sublattice order parameters A^ and Ab, which are k in- 
dependent in the original Hamiltonian, are coupled by C ^ 0, 
and this gives rise to three nonlocal (k-dependent higher par- 
tial waves) order parameters in the transformed basis. Having 
derived the self-consistency equations, next we are ready to 
present main findings of this work mentioned above in the in- 
troduction. 



where we set the Boltzmann constant ks to unity. Fol- 
lowing the usual procedure, we find the lowest-energy state 
of the system by minimizing il with respect to the ampli- 
tudes I A Q I and phase difference $, leading to a set of (three) 
nonlinearly-coupled equations. In addition, we may set the 
total n = n-|- + ni and imbalance p = n\ — nj, number fillings 
using the thermodynamic identities n = —{l/M)dVl/dn and 
p = —(l/M)dn/dh, where < n a < 1. This procedure 
leads to a set of (five) self-consistency equations that needs to 
be solved simultaneously, and four of those can be explicitly 



III. NUMERICAL RESULTS 

In this section, we present our numerical results that are 
obtained by solving Eqs. ©-(O for a self-consistent set of 
|Ayi|, |A_b|, h and \x as a function of given n, g and C val- 
ues. Here, we consider only the ground states of population- 
balanced Fermi gases, and set T = and h = 0. Moti- 
vated by the success of earlier theoretical works on the BCS- 
BEC crossover problem [27-29], we emphasize that the given 



set of four mean-field equations (that are suitably generalized 
here to the checkerboard superlattice model) is expected to de- 
scribe the qualitative physics well for weak interactions at all 
T, and even for moderate and strong interactions at T = 0, as 
long as the single-band tight-binding approximation remains 
valid. In order to isolate the effects of a nonzero C from that 
of checkerboard hopping, let us first analyze the uniform hop- 
ping t\ = t 2 = t case. 
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FIG. 2: (Color online) The superfluid order parameters |Aa| in (a) 
and | As | in (b) (both in units of hopping i) are shown as a function of 
total number filling n for a set of checkerboard potentials C when the 
interaction strength is g — 5t. The particle-hole symmetry implies 
simultaneous transformation of A — > B and \i — s> -/i (or n — ► 
2-n). 



A. Uniform hopping: t\ — ti = t 

In Fig. [2] we set g = 5t and show | A a I an d I A B | as a func- 
tion of n for a set of C values. When G = 0, lA^I = \Ab\ 
is symmetric around half filling (n = 1 or /i = 0), which 
is a consequence of the particle-hole symmetry of the parent 
Hamiltonian, and its maximum value at half filling is a con- 
sequence of the lattice density of states effect. As mentioned 
in Sec. (TTJ when C ^ 0, the particle-hole symmetry implies 
simultaneous transformation of A — >• B and \x — >• — n (or 
n — > 2 — n.). This is clearly illustrated in all of our numerical 
results, and therefore, it is sufficient to restrict our discussion 
only to low (particle) fillings < n < 1. We note that while 



C^ increases lA^I for low fillings as a function of C, it 
decreases |Ab|. This is because since C > lowers (raises) 
the onsite energy of the A (B) sublattice, the local chemical 
potential of the A (B) sublattice is also lowered (raised) by 
C. Therefore, the A (B) sublattice is effectively becoming 
more and more strongly (weakly) interacting as a function of 
C. This effect is clearly seen in Fig. [3} a), where we plot | A^| 
and | As | as a function of C for a set of n values. In the 
C ^> t limit, we expect \Aa\ — g^fn{l — n) and |As| = 
for < n < 1, which is in perfect agreement with our numer- 
ical results. 

In addition to these findings, we find at precisely the half 
filling that the system undergo a superfluid-normal quan- 
tum phase transition beyond a critical C, as illustrated in 
Fig. 0a). To gain intuitive understanding of this transition, 
we analyze the single quasiparticle/hole excitation spectra of 
the interacting system at n = 1, and therefore, set fx = 
and \A A \ = \A B \ = |A| in Eqs. © and ©. This gives 
X — and a doubly degenerate E\^ s — /3 s ^/|ek| 2 + |A| 2 , 
whose form is the same as the C = spectra if we iden- 
tify | A c | 2 = |A| 2 + C 2 . Here, we recall that the single- 
quasiparticle/hole dispersion relation of the noninteracting 
(.9 = 0) system is simply given by eta = 7a \/|ek| 2 + C 2 
(see Sec. Ill All . Therefore, these results suggest that the sys- 
tem does not favor superfluidity and the normal state with 
|A| = becomes the ground state when C > |Ao|, where 
| An | is the C — value. Our numerical results shown in 
Fig-Hb) are consistent with this analysis, where the the criti- 
cal C values exactly coincide with | Ao | along the phase tran- 
sition boundary. Given that |An| = (g/2 — 4t 2 /g) \/n(2 — n) 
in the strong-coupling (g ^> i) limit, by setting n = 1, we ob- 
tain C = g/2 — 4t 2 /g as the asymptotic limit of the boundary, 
and this is in perfect agreement with our numerical results. 
When the critical C -> g/2 > t, we find that the A (B) 
sublattice has ha —> 2 (jib — > 0) so that it corresponds to a 
sublattice band insulator (fully-empty sublattice band) form- 
ing a checkerboard insulator. This intuitive result can also 
be obtained by noting that ua — jib = —(2/M)dQ,/dC = 
Ylks(®Ek. s /dC) tanh[(£'ks/(2T)], which reduces to ua — 
n B = (4/M)J2 k C/^/C 2 + |A| 2 + |e k | 2 at half-filling at 
zero temperature. Therefore, the ground-state of the half- 
filled system first changes from a superfluid to normal when 
C = |A | at which point |A| vanishes, and then the normal 
state rapidly turns into a checkerboard insulator as C/t —> oo. 
For instance, ua — Ub becomes 1.9,1.95,1.99 and 1.999 
when C/t is approximately set to 4.9, 7.4, 17 and 55, respec- 
tively, in the normal state. Since these numbers are indepen- 
dent of g/t, the phase transition is almost (up to one percent 
deviation) directly from the superfluid to a checkerboard insu- 
lator when g/t > 34. 

We note in passing that while having a nonzero C leads to a 
staggered pattern not only in the sublattice order parameters, 
i.e. lA^I 7^ |Ab|, but also in the sublattice number fillings, 
i.e. ua 7^ Ub, we avoid calling the ground-states of the sys- 
tem a supersolid (when |A^| ^ \A B \ and tia ^ n B ) or a 
charge-density- wave insulator (when |A^| = |A^| = and 
riA 7^ n B is an integer number). This is because since C^0 
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FIG. 3: (Color online) (a) The order parameters |Aa| and |Ab| (in 
units of hopping t) are shown as a function of the checkerboard po- 
tential C (in units of t) for a set of total number fillings n when 
g = 5t. (b) The superfluid-normal phase transition boundary is 
shown as a function of the interaction strength g (in units of t) and C 
at half filling (n = 1). Note that C = |Ao| coincides with the tran- 
sition boundary, where |Ao| is the magnitude of the order parameter 
at C = 0. The normal state rapidly turns into a checkerboard insula- 
tor as C/t — > oo, which is not shown in the figure but discussed in 
Sec. lITTAl See Sec. HITBl for the blue-dotted line. 



breaks the translational invariance of the lattice, directly caus- 
ing such an alternating order parameter and filling patterns , 
we believe it is important to distinguish our superfluid and 
checkerboard insulator phases from the true supersolid and 
charge-density-wave insulator ones, for both of which the 
translational invariance is broken spontaneously due for in- 
stance to the presence of nearest-neighbor interactions. 



B. Checkerboard hopping: £i =^ £2 

Before we present our concluding remarks, here we discuss 
the possibility of having an alternating hopping amplitudes. 
Equation (Q]l indicates that the effects of small deviations from 
the uniform hopping, i.e. when \t\ — t-z\ <C t\ + £2, can be 
taken into account (to a very good approximation) via chang- 
ing the normalization factor t that is used in the previous sec- 
tion to (£1 + £2)/2. However, in the asymptotic £2 — ^ limit, 



Eq. ((TJ gives |ek| = 2£i cos(fcid/2), and one may expect 
major quantitative differences. For instance, the superfluid- 
normal phase transition boundary 2C/£i of a half-filled sys- 
tem is shown in Fig. |3jb) as a function of 2g/ti (blue-dotted 
line), and it is clearly shown that the phase boundary devi- 
ates substantially from that of the uniform hopping case. We 
again note that the normal state rapidly turns into a checker- 
board insulator as C/ti — > 00, for which tia — ub becomes 
1.9, 1.95, 1.99 and 1.999 when C jt\ is approximately set to 
4.2, 6.1, 14 and 45, respectively, in the normal state. Having 
discussed the numerical results, next we briefly summarize the 
main findings of this paper. 



IV. CONCLUSIONS 

To summarize, here we studied the ground-state phases of 
Fermi gases loaded into a two-dimensional checkerboard su- 
perlattice potential, i.e. a double-well optical lattice, con- 
sisting of two interpenetrating square sublattices A and B. 
We described this system with a Fermi-Hubbard type lat- 
tice model which includes, in addition to the usual nearest- 
neighbor hopping t and onsite (attractive) density-density in- 
teraction g, a sublattice-dependent local (onsite) energy C. 
Within the single-band tight-binding BCS mean-field theory, 
we reached the following conclusions for such a Hamiltonian. 
First, we showed that the s-wave sublattice order parameters 
A a = \A A \e i<t,A and A B = |A B |e^ B , which are k indepen- 
dent in the original Hamiltonian, are coupled by the presence 
of a checkerboard potential (7/0, and this gives rise to a 
Hamiltonian that has the form of a two-band superfluidity with 
three (two intraband and an interband) nonlocal (k-dependent 
higher partial waves) order parameters in the basis where the 
single-particle Hamiltonian is diagonal. We studied the evolu- 
tion of these order parameters as a function of particle filling, 
interaction strength and checkerboard potential, and found 
that the system always prefers the 0-phase ((f) a — 4>b) solu- 
tions but never the 7r-phase {<f>A = 4>b + i") one. In addition, 
we found at precisely half fillings that the ground-state of the 
system undergo a superfluid-normal quantum phase transition 
beyond a critical C, the threshold of which is precisely deter- 
mined by the magnitude of the order parameter at C = 0, and 
that the normal state rapidly turns into a checkerboard insula- 
tor as C increases. One may extend this work in many ways, 
and motivated by the ongoing experiments Il30l - l34ll . we are 
especially interested in studying the effects of artificial gauge 
fields on the ground-state phases of the system, e.g. the so- 
called optical flux lattices. 
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